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Abstract. We derive a rigorous scaling law for minimizers in a natural version of the 
regularized Cross-Newell model for pattern formation far from threshold. These energy- 
minimizing solutions support defects having the same character as what is seen in exper- 
imental studies of the corresponding physical systems and in numerical simulations of the 
microscopic equations that describe these systems. 



1. Introduction 

This paper reports on some recent progress that has been made in the analytical modeling 
of defect formation, far from threshold, in pattern forming physical systems. We will take a 
moment here to very briefly sketch the physical and mathematical background that motivates 
what is done in this paper. 

The relevant class of pattern-forming physical systems to consider are those in which the 
spatial physical field can be described as planar and the first bifurcation from a homogeneous 
state, having arbitrary translational symmetry in the plane, produces a striped pattern which 
has only a discrete periodic symmetry in one direction. This symmetry-breaking occurs at 
a critical threshold; above this threshold the pattern can deform and, further away, defects 
can form. It is the desire to understand and model this process of defect formation that 
motivates our study. 

A good particular example of these kinds of physical systems is a high Prandtl number 
Rayleigh-Benard convection experiment. The critical threshold in this example is the critical 
Rayleigh number at which fluid convection is initiated from the sub-threshold homogeneous 
conducting state. The " striped pattern" here can be taken to be the horizontal cross-section 
of the temperature field at the vertical midpoint of the experimental cell in which convection 
rolls have formed. 

Because of its periodic structure, the striped pattern can be described in terms of a 
periodic form function of a phase, 9 = k ■ x, where the magnitude of k is the wavenumber 
of the pattern and the orientation of k is perpendicular to the stripes. Here x = (x, y) is a 
physical point in the plane. Even though the striped pattern will deform far from threshold, 
over most of the field (and in particular away from defects) it can be locally approximated 
as a function of a well-defined phase, 9(x), for which a local wavevector can be defined 
as k = V8 which differs little from a constant vector unless one varies over distances on 
the order of many stripes in the pattern. This slowly-varying feature of pattern formation 
far from threshold motivates the introduction of a modulational ansatz in the microscopic 
equations describing these physical systems from which an order parameter equation for 
the behavior of the phase can be formally derived. This was originally done by Cross and 
Newell [3]. These equations are variational and from our perspective it is advantageous to 
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study their solutions by studying the behavior of the minimizers of the variational problem. 
The version of the variational problem that we study corresponds to the following energy 
functional on a given domain Q with specified Dirichlet boundary values. 

(1) ^(0)=^ f (A^Q) 2 dX + - [ (l-\V^e\ 2 ) 2 dX , 

Jn t 1 Jn 

which is expressed in terms of slow variables stemming from the modulational ansatz men- 
tioned above: X = (X, Y) = (fix, fiy) ; = -. 

We refer to this functional as the regularized Cross-Newell (RCN) Energy. It consists of 
two parts: a non-convex functional of the gradient (the CN part) plus a quadratic functional 
of the Hessian matrix of 0, which is the regularizing singular perturbation. Without this 
regularization, the CN variational equations admit non-physical caustic formation. Instead, 
by studying the limit of minimizers of £ M as fi — * 0, one may be able to identify the formation 
of a physical defect as a limiting jump discontinuity or other kind of singularity in the 
wavevector field associated to the /z-indexed family of minimizing phase fields. 

For more details on what has been rather tersely outlined above, we refer the reader to [4] 
where analytical results on the asymptotic limit of minimizers for RCN and their defects in 
certain geometries are also derived. See also [6] where further refinements and generalizations 
are developed. We further mention that the variational problem associated to (1) also arises 
in other physical contexts (unrelated to pattern formation) where it is known as the Aviles- 
Giga energy [2]. 

We now turn to the focus of this paper. The kind of defects that are seen to arise far 
from threshold are not supported by asymptotic minimizers of (1) if the class of functions 
over which one is varying is restricted to be single-valued phases. In particular, one can 
see for purely topological reasons that this restriction rules out disclinations [4]. In [5], 
physical, numerical and experimental arguments are developed which make a strong case in 
support of the hypothesis that the correct order parameter model for the phase in pattern 
forming systems far from threshold should come from a variational problem admitting test 
functions which are multi-valued and in particular two-valued. In physical parlance this is 

— * 

often expressed by saying that the wavefield k should be allowed to be a director field; i.e. 
an unoriented vector field. One figure (see Fig. 1 below) from [5] will help to crystallize the 
issue and the focus of this paper. 

This figure shows seven numerical simulations, each done in a horizontal strip, of a solution 
to the Swift- Hohenberg equation which is a generic model of microscopic equations for a 
pattern forming system. Each of these is run far from theshold but with differing boundary 
conditions imposed at the edges. In each case the boundary conditions impose a constant 
orientation of the stripe at the edges such that the normal to the stripe is (cos(a), sin (a)) 
along the top edge and (cos(a), — sin(ai)) along the bottom edge. The only thing that 
changes from one simulation to the next is the value of a which in the figure is recorded 
on the left in each respective cell. The results of [4] together with symmetry considerations 
establish that for an analogous domain and boundary values, the asymptotic minimizers of 
(1), within the class of single- valued phases, should have the form shown in the bottom- most 
cell of Figure (1). That is, they should have wavevectors very close to (cos(a), sin(a)) in 
the upper region of the cell and very close to (cos(a), — sin(a)) in the lower region of the 
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Figure 1. The "Swift-Hohenberg" zippers. The patterns are determined by 
minimizing the Swift-Hohenberg energy functional for various choices of the 
angle a that detemines the slopes of the stripe patterns as y — > ±00. 

cell with a boundary layer around the mid-line in which the wavevector transitions smoothly 
but rapidly from one state to the other. These minimizers are dubbed knee solutions in 
[4] and in the limit as /i — > 0, they tend to a configuration in which there is a sharp jump 
in the wavevector along the mid-line. This kind of defect is called a grain boundary. In 
other words, the theory for (1) with single- valued phases predicts that the grain boundary 
should be the limiting defect independent of the value of a. The different result appearing in 
Figure (1) was one of the pieces of evidence sited in [5] to argue the necessity for the larger 
variational class of multi-valued phases, even in such simple geometries as those of Figure 
(1). In this paper we are going to carry out a careful analytical study of the RCN variational 
problem in exactly this geometry but within a larger class of two-valued phases. We will 
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firmly establish that the form of the asymptotic minimizers in this more general model does 
in fact depend non-trivially on a. In addition, the construction of test functions in section 3 
and the numerical simulations in section 5 gives some intuitive and experimental support to 
the belief that the stable solutions of the RCN equations qualitatively resemble what is seen 
in the Swift-Hohenberg simulations. In [5], the term Swift-Hohenberg "zippers" was coined 
to refer to the problem studied in Figure (1). In this paper we will be studying Cross-Newell 
zippers. 



2. Setup 

We are given an angle a that determines the boundary conditions on the pattern as 
y — > ±oo by 

V0 — > (cos (a), ± sin(a)) as y — > ±oo. 

Note that this differs from the setup underlying the Swift-Hohenberg zippers in that the 
boundary conditions are placed at ±oo in the y-direction rather than at finite values of y. 
This simplifies our technical considerations in that we don't need to worry about adjusting 
the location of these boundaries as a changes. Also, all of the patterns we want to consider 
here are shift-periodic in the x-direction. This allows us to reduce our study to domains that 
are periodic in x. We introduce the (small) parameter e = cos(a) and we define the period 
I = n/e. We consider the following variational problem on the strip S e = {(x, y)\0 < x < 
l,y>0}: 

Minimize T^\§\ a, 5} given by 

jr[0 ;a ,5] = J J {[A9) 2 + (1 - \V9\ 2 ) 2 }dxdy 

over all a G [0, 1], S G M and 9 satisfying the boundary conditions 
(2) 9(x, 0) = for < x < al; 

9 y (x, 0) = for al < x < I; 
9(x, y) — ex is periodic in x with period I for each y > 0; 

9(x, y) - \ex + y/l - e 2 y + s\ G H 2 (S e ). 

We take a moment here to explain the considerations that have motivated the mixed Dirichlet- 
Neumann boundary conditions here, the first two boundary conditions in (2) above. We 
argued in the introduction that in order to capture the physically relevant minimizers, the 
RCN variational problem needed to allow for multi-valued phases in its admissible class of 
test functions. However, the numerical results on the Swift-Hohenberg zippers suggest that 
in certain symmetrical geometries the appropriate multi-valuedness can be introduced in a 
tractable fashion. Indeed in the case of the SH zippers we see that the symmetry of the 
boundary conditions between the upper and lower edges of the domain is preserved in the 
symmetry of all of the exhibited solutions about the middle horizontal axis; i.e., the reflection 
in y about the y = axis. This suggests that a single-valued phase could describe the solu- 
tion in the upper half-plane with the solution in the lower half-plane given as a symmetric 
reflection of that in the upper half-plane about y = 0. 

Figure 2 illustrates two instances of the form that we expect these zippers to take in the 
infinite (in y) geometry. The figure on the left illustrates level curves {stripes in the parlance 
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Figure 2. An illustration of the appropriate boundary conditions. 

of the introduction) of what we will shortly define to be a self-dual knee solution. This is 
indeed symmetric about the mid-axis, which we will take to be the y = axis; moreover, 
one can see that its gradient field along y = is tangential to this axis. Thus the gradient 
field in the upper half-plane is completely symmetrical to that in the lower half-plane under 
reflection about y = 0. 

However, for the striped pattern on the right in figure 2, this is not the case. There are 
regions, illustrated for example by the darkened interval along y = 0, where the gradient 
field is tangential to this axis; but, there are other regions, illustrated for example by the 
lightened interval along y = 0, where the gradient field needs to be perpendicular to this axis. 
By reflection symmetry this field will point upwards in the upper half-plane and downward 
in the lower half-plane. This cannot be supported by a vector field but it is allowable for a 
director field. This indicates that in this region a two-valued phase is required. 

To get at the conditions on the phase itself we observe that patterns of the type illustrated 
here are analytically given in terms of a form function F of the phase 9 = 9(x, y) such that F 
is locally periodic of period 2tt in 9 and such that F(9(x, y)) is even in y and smooth in (x, y). 
In order to allow 9 to be two-valued we also require F to be even in 9. (An example of a 
global form function having these properties is F = cos.) It follows from these requirements 
that either 9(x,y) is even in y, in which case 9 y (x, 0) = 0, a Neumann boundary condition; 
or, 9(x,y) is an odd function of y modulo n, in which case 9 y (x,0) = nn for some integer 
it, a Dirichlet boundary condition. Thus to realize the pattern on the right in figure 2 in 
terms of a single-valued phase in the upper half-plane, we would need to take the Neumann 
boundary condition on the darkened interval and the Dirichlet boundary condition on the 
lightened interval. This is what we have done in (2). For the self-dual knee pattern on the 
left we would take the entire boundary condition to be Neumann. 

The functional J- e is the RCN energy functional but with the scaling /i removed. It is 
appropriate to do this because the demonstration that the nature of the RCN minimizers 
depends on a is independent of this scaling. The first and the second conditions impose 
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a mixed Dirichlet/Neumann boundary condition at y — 0, the third condition imposes 
(shifted-) periodicity in x and the last condition ensures that the test functions 9 approach 
the straight parallel roll patterns ex + — e 2 y + 5 as y — > oo. 

Note that the dependence of the functional T e on the parameter e is through the depen- 
dence of the domain S e and the boundary conditions on e. The parameters a and 5 are 
determined by minimization. The parameter a is a measure of the fraction of the bound- 
ary at y = that has a Dirichlet boundary condition, and 5 represents the asymptotic 
phase shift, that is the difference in phases between the test function 9 and the roll pattern 
0(x, y) = ex + VI - e 2 y which satisfies 0(0, 0) = 9(0, 0) = 0. 

The case where a is set to zero is considered in earlier references [4]. The test functions 
9(x, y) satisfy a pure Neumann boundary condition at y = and the minimizers in this case 
are the self-dual knee solutions 

9 n eu(x,y) = ex + log(cosh(Vl - e 2 y)). 

These solutions have an asymptotic phase shift of — log(2) and the energy of the minimizers 
in the strip S e is given by 

(3) ^[9 neu ,0,-\og(2)] = AnV ^ . 

The existence of (9 e , a e , S e ) minimizing jF e can be shown from the direct method in the 
calculus of variations. We also prove the following results about the minimizers, and their 
energy - 

Theorem 2.1. Upper bound 

There is a constant E such that J re [9 t ; a e , 5 e ] < E for all e G (0, 1]. 

We prove this result in sec. 3 by exhibiting an explicit test function satisfying this bound. 
Note the implication that the minimizers for sufficiently small e cannot be the self-dual 
solutions, since the energy in Eq. (3) diverges as e — > 0. Consequently, a e > for sufficiently 
small e. 

Theorem 2.2. Lower bound 

There are constants E\ > and e > such that, even for the optimal test function 9 e and 
the optimal parameter values a e and 5 e , we have J re [9 t ;a t ,S e ] > E\ for all e < e . Further, 
there are constants < ct\ < a<i such that 1 — a^e < a e < 1 — a±e for sufficiently small e. 

We prove this result in sec. 4. Combining this result with the preceding theorem, we 
obtain a rigorous scaling law for the energy of the minimizer, and for the quantity (1 — a) 
as e — > 0. As a corollary to Theorem 2.2, we find that an 0(1) part of the energy of the 
minimizer concentrates on the set, al < x < I, < y < 1. This can be interpreted as saying 
that a nontrivial part of the energy of the minimizer lives in the region of the convex- concave 
disclination pair [5]. 

3. Upper bound 

We will first show an upper bound for the enrgy functional uniform in e, by constructing 
a family of explicit test function whose energy is uniformly bounded. The idea for the 
construction of these test functions comes from the self-dual ansatz [4] which requires that 
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the energy density of the functional T should be equi-partitioned between its two terms. 
Functions satisfying this ansatz solve the self-dual (resp., anti-self-dual) equation: 

(4) A6 = ±(1-\V6\ 2 ). 

Solutions of this equation can be constructed via the logarithmic transform 

9 = ± log u 

which reduces (4) to the linear Helmholtz equation (5). We refer the reader to [4, 6] for more 
background on self-dual reduction. 

3.1. Self-dual test functions for the CN-Zipper problem. 

3.1.1. Existence. We consider the Helmholtz equation in the upper half-plane, 



(5) Am - u = 
subject to the mixed boundary conditions 

(6) u(x,0) = e _mr n£ < x < (n + a)£ 

(7) u y {x,0) = {n + a)£ < x < (n + 1)1 



and with asymptotic behavior for large y given by const. exp(— ex — Vl — e 2 y) where I = ir/e 
and a G (0, 1). 

We seek a shift-periodic solution, meaning that we change variables tow = e ex u(x,y) and 
look for periodic solutions of 

(8) Lw = Aw - 2ed x w - (1 - e 2 )w = 0, 

with boundary conditions of periodicity in x of period £\ mixed boundary conditions at y — 0, 

(9) w(x, 0) = e ex 0<x<a£ 

(10) w y {x,0) = a£<x<£; 

and with asymptotic behavior for large y given by const. exp(— \/l — e 2 y). Given such a u, 
9 = — logw would satisfy the boundary conditions (2). (However, for notational simplicity, 
in the remainder of this section we will set 9 = logw.) 

We now let S € denote the half-cylindrical domain, ^-periodic in x and with y > 0. The 
existence of a weak solution to (8) satisfying the above boundary conditions can be estab- 
lished via the Lax-Milgram theorem with appropriate energy estimates. However, in order 
to derive uniform asymptotic energy estimates (as e — > 0) for the CN Zipper problem we 
need to go beyond existence results and try to construct a more explicit representation of 
the solution to (5). Unfortunately, at present, the solutions one can construct using Greens 
function methods and the like do not yield sufficient a priori boundary regularity near y = 
to control the asymptotic behavior of the energy in this finite part of S € . We will there- 
fore instead study solutions of a self-dual problem with modified boundary conditions (more 
precisely, with pure Dirichlet boundary conditions). Subsequently we will make a local mod- 
ification of these solutions near the boundary to produce functions (no longer global self-dual 
solutions) whose asymptotic energy we can control and which are valid test functions for the 
Cross-Newell Zipper problem. 

The modified boundary value problem we consider is (8) with (9-10) replaced by the pure 
Dirichlet boundary condition 
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(9') . m _ J e ex < x < at 

(10') W[X ^ } ~ \ q a (x) a£<x<£ 

where q a (x) is a function which smoothly interpolates, up through second derivatives, be- 
tween e ex at x = a£ on the left and e ex ~ n at x = £ on the right. There are clearly many 
choices for such a function; the precise choice for our purposes will be made later at the 
end of subsection 3.2. By elliptic regularity [7], the solution to this boundary value problem 
satisfies w(x,y) G H 2 (S e ). In the following sections we will construct the solutions to this 
problem and study its asymptotics relative to the RCN energy T e . 

3.1.2. Explicit Construction. The whole plane Green's function for the Helmholtz equation 
(5) is explicitly given in terms of the Bessel potential [7] : 

(11) G(x, y; t, v ) = ± J™ e^j exp (-!{(* - if + (y - ??) 2 }) . 

In terms of this Green's function we can then represent a solution to (5), with asymptotic 
behavior for large y given by const. exp(— ex — yl — e 2 \y\), as 

/oo 
p e (0G(x,y; £,0)de. 
-oo 

Note that 

/oo 
-oo 

/oo 
ff(t)G v (x,y;Z,0)dt 
-oo 

solves (5) with respect to the standard Dirichlet boundary condition which equals minus the 
jump of u € y along the x-axis. One may check directly (see (22)) that in fact p e (£) = — 2m^(^, 0) 
almost everywhere. Integrating (13) with respect to y gives 

/oo 
p%t)G{x,y-^)dt. 
-oo 

Since both u e (x,y) and the RHS of (14) decay as y j oo, it follows that f(x) = 0. This is 
consistent with the ansatz (12), taking p e (£) to be the jump in the normal derivative of if 
along y — 0. 

We make the following shift-periodic ansatz for p e , 

p *(t + e)e«t+Q = p<(0e*. 
With this one can expand out (12) more explicitly as 
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u e (x,y) = 

(17 ) - ^l^^V«(flZe-p(^). 

(18) = * g -«^»& /W(^x>p ( «? - + »^» 2 ) 

In (15) we have interchanged the order of integration which is justified by Tonelli's Theorem; 
in (16) we've made the substitution £ = + n£ and in (17) we've commuted the sum past 
the integrals which is justified by monotone convergence-all terms in the series are positive 
and hence the partial sums are monotonic. In (18) we write each summand as a single 
exponential and then appropriately complete the square in each exponent. Finally in (19) 
we apply Jacobi's identity [9]. Here $3 is one of the Jacobi theta functions, in this setting 
explicitly given as 

/on , q f-x + 2et -4ttA i o ^ -faoVt (2ixn ( 2nt\\ 

(20) rf^—^—.—j = l + 2ge(,)-cos^- r ^--jj 

Finally, from (19) we can express our candidate for the solution to (8), (9' — 10') as 
where p e (£) = p e (£)e e ^. 

3.1.3. Data Characterization, periodized and in Fourier Space. From the previous sections 
we have that p e (£) is periodic of period l\ also w e (x, y) is periodic in x of period £ and = e ex 
along (0, a£) when y — 0. 

Moreover, taking the Fourier transform of (21) one finds that the Fourier coefficients, in 
x, must satisfy 

(22) {«^)}(n,y) = 1 1 .. 2 M0}We-V^(^ 

2 a/1 + e 2 (2n + ?) 2 

for each value of y. Taking the limit as y — > on both sides of (22) gives 

(23) = 5 7l , 2 ! , -v lMOK")- 

^ a/1 + e^(2n + %y 

This is a determining conditions for p e (£). We note that differentiating (22) with respect to 
y and setting y = demonstrates that p e (x) = —2Wy(x, 0), at least in the L 2 sense. 



Since w e (x,0) G H 2 ^ 1 ), it follows, by comparison, that 2y/l + e 2 (2n + i) 2 w e (n) G h l (Z). 
Given this we can now define 

(24) p e (x) = {2^/1 + e 2 (2n + i) 2 ^(n)} V (x) 

which characterizes p e as an element of H 1 (S 1 ). It follows from Sobolev's lemma [7] that 
p e can be taken to be continuous. This last observation also justifies the existence of the 
Fourier coefficients {p e }(n) that were formally introduced in (22). 

3.1.4. Large y asymptotics. We now determine the large y asymptotics of (12). By (22), w e 
has a Fourier representation given by 

= w^)e-^y + O (V^w) . 

(We note that for large y this series converges uniformly to a smooth, in fact real-analytic, 
function of x.) Moreover, w e (0) = | J* w e (x, 0)dx is non-zero since by the maximum principle 
[7] applied to the elliptic PDE (8) on the cylinder [0, £] x (—00, 00), the integrand, w e (x, 0), 
is non-negative and in fact, by (9' — 10'), non-vanishing on [0, £] (the definition of q a which 
we give later will insure that this is so). 

3.2. Energy Estimates. We will now try to show that the regularized Cross-Newell energy 
of 9(x,y) = \ogu(x,y) is uniformly bounded in e. This would establish a uniform (in e) 
upper bound for the energy minimizers. Recall that the energy is calculated by integrating 
the energy density over the domain S e . Making this estimate breaks naturally into the 
consideration of two regions: [0,£] x {y > M e } and [0, £] x {y < M t } where M e is to be 
determined. 

We remark that the so-called "knee solution" of the self-dual equation provides an upper 
bound for the energy for values of e bounded away from zero. So we only need to be concerned 
with small values of e. Since u e (x,y) = e~ ex w e (x,y) solves the Helmholtz equation, it will 
suffice to bound the density (1 — \V9 e \ 2 ) 2 (since the integral of this density equals that of 
(A6 e ) 2 for self-dual solutions). 

3.2.1. Estimates in [0, £] x {y > M e }. We begin by considering the domain for large y. Since 

Ve e (x,y) = — r (x,y)=( \+—-( x ,y), 

we may reduce our considerations to studying the asymptotics of w e and its first derivatives. 
It will be convenient to replace the convolution integral in (21) by the Fourier series whose 
coefficients are the product of the Fourier coefficients of p e and the $3 series. This results in 
the following alternative representation of w e : 

w e {x,y) = -j=^ / — e V 'VWe < • 
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With the change of variables, 

t 

s — - 

y 

this representation takes the form 

(25) w (x,y) = J — 2_^ —e V s " J p\n)e 



47T 



2-Kinx 

I 



s 2 



where s n = 2 ^ 1+ ^ 2 ^ e cr ^^ ca ^ P om t °f the exponent is s = s n and the expansion of 
the exponent in the nth term of the series near this critical point has the form 

S + l = 2_( 1+ (s-Sn)* +o r(8- 8n » 



An asymptotic expansion in large y may be developed for the integral in each term of the 
series (25) by the method of Laplace. By the uniform convergence of the series (for large y), 
the asymptotic expansion of the series is equivalent to the sum of the asymptotic expansions 
from each term. We implement this strategy to find the leading order, large y behavior, and 
next corrections, for w e , w e x and w y : 



w'(x,y) = JH-Yde-Z r^^e-^(^ 3 ^(n)e^, 

, , . _ . / y V \ i V f°° dZ V I J2 i /nr~3\\ -~ , . 2winx 

wi(x,y) = -2ie* -f-y sle * / r e ^ 2 +U(Z >)np e (n)e~ ~, 

, , v 1,1/?/ \ — i _JL f°° dZ V—( 'A-L.n!'/ A \\ ^~ , \ 2-Kinx 

w e y (x,y) = — » £ -^ t fV% 2 e «» / + ^^) p e( n ) e — , 

ly I V 47r ~Tn, J-i (l + z)2 



where in the n th term of each series, z = ^zSji respectively. We can now apply Laplace's 
method to each term and then observe that the dominant contributions for large y come 
from the 0, +1, —1 Fourier modes. Retaining just these we derive the following asymptotic 
behavior for Vlogw e : 



-±(x,y) = ^^e-^(P(l )e ^)^(ee-^) 



^x,y) = - v'F^ Vp(0) ' K 1 



= + ±- + O (e 2 e- 2 ^) . 

2y V J 

Based on these asymptotics we can now estimate the energy in the large y domain. 



W = 



-x/r^2 



O ( J + e 2 e- 2 ^) 
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from which it follows that 



|W| 2 = l + 0(eV^)+0^ 



and so the energy density 



(l-IWI 2 ) 2 = ( £ V-) +0 ge-^)+0(i). 

^From this it follows that the "large y" part of the total energy is bounded as 

T € < 

Thus, if we take M e = c/e, this part of the total energy will remain finite as e — > 0. 

3.2.2. Estimates in [0,£] x {y < M e }. We next turn to consideration of the energy density 
in the finite part of the domain. To facilitate this consideration we will sometimes make the 
uniformizing change of variables z = e£ and h = ex in the Jacobi theta function (20): 

(x-£) + 2et -4irt\_ Q f-{h-z) + 2e 2 t -4eV 
1 



(26) *,P^.^ =4 



7T 7T 



In what follows we will assume that a is chosen to depend on e in such a way that 1 — a e = 
0(e). 

We will make use here of the following single-layer potential counterpart of the double- 
layer potential representation (21), which in fact can be deduced directly from a change of 
variables in (25): 

w e (x,y) = 



We study the asymptotic behavior of the convolution integral in (27) for x G (0, at) and 
for times t of order less than 1/e: 



7T Jo V T T / 

1 r/ z \ f-{h-z) + 2eH -4e 2 A 

= e ex + o(e), 



where in the third line above, the form of the integrals follows from making the change of 
variables as in (26). In the second integral we smoothly extend q a {~) — & z to be zero on 
(0, air). The final line follows for t of order less than 1/e because in this regime the Jacobi 
theta function inside the convolution behaves as a Dirac comb as e — > 0. The second term has 
this asymptotic behavior because h € (0, air) and the support of ? (f) — e z is complementary 
to this interval, so that this integral decays exponentially to zero with e, as with a Dirac 
sequence away form its support. 

Based on (28) we can estimate w e as 



w e (h,y) = 



ti 



-((l-e 2 )t+^) (e e 



+ o(e)) 



(29) 



= e ex e- VT ^ y + o(e). 



The evaluation of the previous integral may be deduced from a basic Bessel identity (see [1] 
9.6.23). 

In order to estimate V8 e , we also need to estimate the x and y derivatives of w e (x, y). To 
this end we first consider the x-derivative of the internal convolution integral which equals 



(30) 
(31) 



1 



d x - J Q ( qa (- e ) - e ) ^ { , | dz 



-c- I" U(-)-e z ) d z d 
7T Jo v e J 



7T 

-(h -z) + 2eH -4eH 



dz 



IT IT 

-{h-z)+2e 2 t -4e 2 t\ . 
3 i , I dz. 



IT 



7T 



Integrating by parts, the above derivative may be rewritten as 



(32) 



- / ed z (q a (-) ~ e 2 
7T Jo V e 

o(e) for h G (0, air), 



-(h-z) + 2e 2 t -AtH 



IT 



IT 



dz 



as for the second integral in the last line of (28). Thus, 



u c 



(33) 



-e + 



-e + 



OO rft ((1 _ e 2 )t+ ^) 1 



v 7 ^ JO 



o(e) 



e ex + o(e) 



0(e) 



For the y logarithmic derivative we have 
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1 




y 



1 



/o °^- ((1 - e2)t+ « ) l/o^O)^( 

l/ oo fe-« 1 -^)(e- + O (6)) 



(g-Q+get -47rt 




* e -((i-* 2 )t+£) ( e - + o(e)) 



(34) 



1 



The last equivalence follows from a Bessel recurrence identity [1], formula 9.6.26, together 
with formula 9.6.6. 
Thus we finally have 



and hence (1 — |V6> € | 2 ) 2 = 0(e 2 ). Since the domain [0,a£] x {y < M e } has dimensions 
1/e x 1/e, the total energy in this region is also asymptotically finite. 

3.2.3. Modification of the Self-dual Test Funciton. It remains to estimate the energy in the 
region [a£,£] x {y < M t } which has dimensions 0(1) x 1/e. The question of the finiteness of 
the energy of the 6 e we have been considering in this region is beside the point for general 
purpose of but, this self-dual solution does not satisfy the boundary condition (10) in this 
region. 

As stated earlier we are going to modify the self-dual test function in this region so that 
the boundary condition (10) is satisfied. To that end we fix a small value of 5 and let B(5) 
denote the ^-neighborhood of [a£,£] in S £ . We modify w e in this neighborhood as follows. 
Define 

(36) w £ (x,y) = </) 1 (x,y)w e (x,y) + (f) 2 (x,y)w 2 (x,y) 

where {0i, 02} is a partition of unity subordinate to the cover of S e given by 



and w 2 (x, y) = w e (x, 0) cosh (y), where w e (x, 0) here is defined as in (9' — 10'). One has 



(35) 




Z7i = S e \B(5/2) 
U 2 = B(6) 




1 C/i\B(5) 
B(5/2) 




1 B(S/2) 
C/i\B(5) 



and 0i + 02 = 1. 
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It is straightforward to check that w e (x,y) satisfies the boundary conditions (9) and (10): 



\imw e (x, y) = <fi\{x, 0)w e (x, 0) + <j)i(x, 0)w e (x, 0) 

= (0!(a;,O) + 2 (x,O))w e (x,O) 
= w e (a;,0) 
= e tx 

for x G [0, a£\. 
For x G [a£,£], 

\imw e y (x, y) = ((/>i y (x, 0) + 4> 2y (x, 0)) w e (x, 0) + (f) 2 (x, 0)w e (x, 0) sinh(0) 

y— >0 

= (0i + 2 ) y (x, 0)w € (x, 0) + 
= 0. 

Thus, log w e is an admissible test function for the regularized Cross- Newell variational prob- 
lem. We can now estimate the energy of this test function in B(<5). The energy density in 
this region is bounded and therefore the energy in B(<5) is finite. 

3.2.4. Estimates for the "outer" solution in ([a£,£] x {y < M e }). It remains to estimate the 
energy in {[a£, i\ x {y < M e })\B(S). To proceed with this we will need a more specific 
definition of q a which we now give. 

Note first that by our assumption that 1 — a e = 0(e), the interval [a£,£] remains of size 
(9(1) for arbitrarily small values of e. We will now further pin this down by setting 1 — a e = ce 
for a value of c that is fixed, independent of e. Consequently, [a£, £} is now an interval of 
fixed length en which can therefore also be represented as [£ — en, £]. Recall that q a needs to 
be built so that on this interval it matches, through second order, to e ex at the left endpoint 
and similarly to e^" 71 " at the right endpoint. Toward this end we observe that the required 
leading order value on the right is 1, independent of e while on the left the leading order 
value limits to the stable value of e w as e — > 0. 

Choosing a value v > that is small with respect to C7r, we define a compressed tanh- 
profile that interpolates between the point (xo,yo) = [£ — en + u, e n + 7) and the point 
(xi,yi) = (£ — u, 1 — 7) and where 7 > is another chosen value required to be smaller than 
1. (This last requirement will insure that the positivity claim made at the end of subsection 
3.1 holds.) Explicitly this tanh-profile is given by 



(X - (£ - V)) (X - (£ - C7T + v)) 

Note that the profile of T(x) is independent of e. The only way in which T depends on e is 
that this profile translates uniformly with £ as e changes. We will define q a (x) = T(x) on 
the subinterval [xq, x±] = [£ — cir + u, £ — v\ of [£ — en, £}. 

Next we will define the piece of q a (x) on the left that interpolates between the point 
(£ — C7r, g 71 " - eC7r ) a nd the point (x , y ). Choose a value a > that is small with respect to 
v. Consider the covering of [£ — en, £ — en + v\ by the two sets V\ — [£ — cn, x — a) and 
V 2 — (£ — en + a, x ] and let {ipi(x), ^{x)} be a partition of unity subordinate to this cover 
which means, in particular, that 
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i>2 



1 [£ — C7T, £ — en + a) 

(x - a, x ] 

1 (x - cr, x ] 

- C7T,£ - C7T + (J). 




^ 1 (x)e^ + ^ 2 (x)(e 7r + 7 ). 



It is straightforward to see that with these choices q a {x) is smooth throughout [£ — civ, t — v) 
and satisfies the smooth matching conditions on the left. Moreover, it is clear from the 
functions comprising (37) that q a and its derivatives remain bounded on [£ — en, £ — v) as 
e — > 0. A similar construction may be made on the right; i.e., on {£ — ax + v,£]- This 
completes our description of q a {x). 

The study of the convolution integral in (27) in the region where x G [a£, £} now proceeds 
similarly to what was done in (28) and subsequent formulae. In particular, the analogous 
result to (28) is that for x G {a£, £) and for times t of order less than 1/e: 



with q a (x) here bounded away from zero, independent of e, by our earlier choice of 7. Hence 
the denominators in the estimates analogous to (33) and (34) are under control. In the 
subsequent formulae the roles of e ex and q a {x) are effectively interchanged as above and all 
proceeds as before. The result is that the energy in ([a£, £] x {y < M t }) \B(5) is asymptoti- 
cally bounded like 0(e). 

It thus follows that the total energy of our family of test functions is uniformly bounded 
in e. 



Following the ideas of Jin and Kohn [8], we will prove ansatz-free lower bounds for the 
functional T e by identifying vector fields ^(V^) such that 



(38) 




q a (x) + o(e) 



4. Lower bound 
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This allows us to obtain information about the energy T e purely in terms of the boundary 
conditions on 9. To avoid the proliferation of symbols, here and henceforth, C,C',Ci, etc 
denote (finite) constants whose precise value is unimportant, and different occurrences of 
the same symbol might denote different values of the constants. e±, e 2 , K, K ± , etc denote 
constants that have the same value in all their occurrences. 

Definition 1. A smooth vector function E(p, q) = (Ei,E 2 ) is subordinate to the energy if 



(39) 



9Si(p,g) 



dp 



+ 



9Ei(p,g) + <9E 2 (p,g) 



dq 



dp 



+ 



dq 



<C\\-p 2 -q 2 



for some C < oo. 

If E is subordinate to the energy, it follows that 



|V-E(V0)| < 



dp 



\9xx\ 



dZifaq) + 9E 2 (p,g) 



dq 



dp 



9xy I 



dZ 2 (p,q) 



dq 



'yy\ 



< C\l-6l-6 2 y \\VV6\, 
where we use the identification p = 9 x ,q = 9 y and | VW| 2 = 9 2 xx + 29 2 y + 9 yy . Consequently, 

F e [0;a,8\= [[ {[VW] 2 + (1 - \V9\ 2 ) 2 }dxdy-2 [ 9 x d9 y 

J JS^ ' JdS* 



>2 // |VV0||1- \W9\ 2 \dxdy-2 
2 ' 



9 x d9 y 



> 



C 



[[ V • Y,{V9)dxdy - 2 I 9 x d9 y 
J J V • Y,{V9)dxdy . 



In obtaining the last equation, we use the fact that 



9 x d9 y — 



for the boundary conditions in (2). 

Lemma 4.1. There are constants e 1: K x > such that for all e e (0, l],a e [0, 1] and 5gM, 
we have 



T'\9- a, 5] > 



I -a) 



K x e 2 . 



Proof. Let > be a smooth, compactly supported function such that 

0(0) = 1, 
0(1) < 0(0), 

f(p) — P4>(p 2 ) has a single maximum at p — 1. 
An explicit example of a function whith these properties is 



0(p) 



exp 




1 1 v 

2 (2-p)(p+l) 4 
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pe (-1,2) 

otherwise 



Let b — (1 — a)/e. Define the vector field E(p, q) by 
E 2 (p,( Z )=p0(6 2 p a ) 



Si(p,g) 



[0(6 2 (1 - ?? 2 )) + 2(6 2 (1 - 77 2 )0'(& 2 (1 - ?? 2 ))] ^. 



Since has compact support, it follows that S is bounded on R 2 . An explicit calculation 
shows that the quantities £i iP and £ 2;? are zero. Also, 



d^( P ,q) + dE 2 (p,q) 



dq 



dp 



= |0(6V) + 2&W(&V) - <K^ 2 (i - <? 2 )) 

-26 2 (l-gV(6 2 (l-g 2 )| 



where 



C = sup 



x,y 



<Cb 2 \(l-p 2 -q 2 
(x) + 2x<j)'{x) - 0(y) - 2y<t>'{y) 



x-y 



< sup 1 3(f)' (z) + 2z(f)"(z) 



is clearly finite since is compactly supported and twice differentiable. This proves that E 
is subordinate to the energy. 

Since V • £(V0) = (£2,35 + £i, g )^x» we obtain 



^ V • E(V0) dxdy < C6 2 ^ 

C6 2 



< 



(40) 



2 

C6 2 



l-9 2 x -9 I y )9 xy \dxdy 
J J (I -61- 9 2 y ) 2 dxdy + JJ [VV9] 2 dxdy 
F e [9;a,5] 



Integrating by parts, we have 

J J V • E(V0) dxdy =£ 2 (e, Vl -e 2 )^ - J ' S 2 (0, ^(x, 0))d. 



7r/e 



£ 2 (^(x, 0),0)dx, 

ct7r/e 

where the contributions from the boundaries at x = and x = ir/e cancel due to the 
periodicity. By construction, S 2 (p, 0) = and S 2 (p, q) has a maximum value at p = 1/6. 
Consequently, 



I 



E 2 (0,9 y (x,0))dx = 



r w .o),oM,<^^ 



7T0(1). 



Also, 0(0) = 1 and is Lipschitz so that 



S 2 (e, v 7 !^ 2 ) = e0(6 2 e 2 ) = e0((l - a) 2 ) > e(l - C"(l - a) 2 ), 
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for some finite C Combining these estimates with (40), we obtain 

2tt 



(41) 



F e [0; a, 5} > 



Cb 2 



0(1) _ c"(l - a) 2 ] , 



and rewriting b in terms of a and e yields the desired conclusion. 



□ 



The above lemma shows that the energy grows without bound as the quantity (1 — a) 
becomes small. However, we do not have a priori control on the size of (1 — a). Consequently, 
to obtain a lower bound for the energy, we need a complementary estimate which shows that 
the energy grows as the quantity (1 — a) becomes large. 

To prove this result, we first construct a vector field S subordinate to the energy functional 
as follows: 

Let ip > be a smooth, compactly supported function such that 

V>(o) = i 

PCX) 

/ (i-mm=o 

Jo 

We can always construct such a function, given x > 0, a compactly supported function 
with x(0) = 1. Observe that 



where A , A 1 > 0. Consequently, by an appropriate choice of 77, we get ip(x) = x( x /v 2 ) with 
the required properties. 

We define the functions ((q 2 ) and cr(q 2 ) by 

cq 

(<? ~ rj)ip{rj 2 )drj 



(42) 



C(s 2 ) 



(q-r))(l - r] 2 )^(r] 2 )dr] 



Note that the functions ( and a are well defined for positive values of their arguments, that 
is the expressions on the right hand sides of the above equations are even functions of q. 
From these expressions, we have 

d 2 



dq 



2 C(<7 2 )=#? 2 ); C(o) = o 
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dq- 



<q' 



;i-?W); cr(0) = 



We will also use the same letters ( and a to denote smooth extensions of the functions 
defined above to all of R. We will pick extensions such that the supports of ( and a are 
contained in [— l,oo). 

Let if > be a compactly supported function and set 



(43) 



V(p,q) = <p(p 2 ) [a(q 2 ) - p 2 ((q 2 )] 



2 (£V(£ 2 )) 
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V is now an even function of p and q. Define the vector field £ by 
(44) 

From (44) it follows that 



Z(p,q) = [- ir V,—V 



dp dq 
dHi(p,q) <9S 2 (p,g 



dq 



+ 



dp 



0. 



Also, 



<9£ 2 (p,g) 




d 2 V(p,q) 


dq 




dq 2 



W{p 2 )^{q 2 )\\l-p 2 -q 2 \<C 1 \l-p 2 -q 7 



With C\ = sup \tp(p 2 )ip(q 2 )\ < oo. Finally, an explicit calculation yields 



dXi(p,q) 



dp 



(45) 



dp 2 
< C 2 \l-p 2 



(a(q 2 ) - (7(1 - p 2 )) 



P 2 ip(p 2 ) 



(C(g 2 )-C(i-P 2 )) 



where C 2 can be bounded in terms of the support of f and the maximum values of \<p\, \ <//], \<f"\, \('\ 
and \a'\. Clearly <p and all it's derivatives are uniformly bounded since it is smooth and com- 
pactly supported. From (42), we have 



ip(r] 2 )dr] 



a'(q 2 ) = ±- f '(1 - ITW)*7 



Since ^ is compactly supported, these derivatives vanish as q 2 — > oo, implying that (' and 
a' are bounded for all positive values of the argument. Since a and ( are smooth and are 
identically zero if their arguments are sufficiently negative, it follows that C 2 < oo. It thus 
follows that the vector field £ is subordinate to the energy functional. 

For future use, let us record a few observations that follow directly from the construction: 

Observation 1. S 2 (p, 0) = since S 2 is an odd function of q. 

Observation 2. 

£ 2 , 9 (0,g) = ^(0,g) = V^ 2 )(l-g 2 ). 
Consequently, the non-degenerate critical points are at q — ±1. Differentiating in q, we get 

so that S 2 (0, q) has a maximum at q — 1 and a minimum at g = — 1. 
Finally, 

S 2 (0, g) = 2qa'(q 2 ) = f (1 " ^WK, 
so that S 2 (0, 1) > and S 2 (0, g) — > as q — > oo. 

M = /^(l — £ 2 )ip(€ 2 )d£ will denote the maximum value of E 2 (0, g). 
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Observation 3. 



S 2 (e, vT^) = y?(e 2 ) /^(l - e 2 - £>(£ 2 K 

./o 

> M-Ke 2 

for a constant K < oo. In obtaining the last line, we use 

\1 - cp(e 2 )\ < C.e 2 



JVT~e 

for some bounded constants Ci, C 2 , C 3 



<K£ 2 R 



<^2 
<C 3 6 2 



Lemma 4.2. There are constants e 2 , if 2 > such that for all e G (0, 1], a G [0, 1] and 5 G M, 
we have 

^[6;a,5]> e -^^-K 2 e. 



Remark. For the case a = 0, corresponding to the self-dual minimizers, this estimate 
captures the right scaling of the minimum energy as e — > 0. 

Proof. The proof of the lemma follows from estimating a lower bound for the functional 
T e \Q\ a, 5} using the vector field S that we constructed above. 
For the vector field £ we have 

J J V • £(V0) dxdy =S 2 (e, Vl-e 2 )- - y S 2 (0, 9 y (x, 0))dx 

- / 7 Z 2 (9 x (x,0),0)dx. 

J aw/t 

As before, the contributions from the boundaries at x = and x = nfe cancel due to the 
periodicity. By construction, S 2 (p, 0) = and S 2 (0, q) has a maximum value M at g = 1. 
Consequently, 

/ £ 2 (O,0j,(x,O))da; < 
Jo e 
From observation 3, we obtain 



Since £ is subordinate to the energy, 
(46) 

which yields the desired conclusion. 



F £ [0;a,5] > — [l-a-Ke 2 ] , 



□ 



We can now prove theorem 2.2 using lemma 4.1 and lemma 4.2. 
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Proof. Let b denote the quantity (1 — a)/e. From lemma 4.1 we get 

F € [9;a,6] > ^-Kxe 2 > 3 (e 1 e 2 2 ) 1/3 - 2e 2 b - K x e 2 

where the last inequality comes from linearizing the convex function e\b~ 2 at b = {ei/e 2 ) 1 ^. 
Combining this estimate with the conclusion of lemma 4.2, we get 

F\Q- a, 5} > max (3 M) 1/3 - 2e 2 b - K.e 2 , e 2 b - K 2 e) > (e^) 1 / 3 - 2K ^ e + K ^\ 
If we set 

. ( (e ie 2 y/s ( ei e 2 )^ 

e* = mm ■ , 

V v 7 ^ 2K 2 

for all e < e*, all a G [0, 1] and all 9 satisfying the boundary conditions in (2), we have 

Combining the upper bound J re [9 t ; a e , 5 e ] < E in theorem 2.1 with the lower bounds for 
JF £ in lemma 4.1 and lemma 4.2, it follows that for 

e < min 

we have 

ei 1 - a e 2E 
< < 




2E t e 2 
Consequently, 

1 — a 2 e < a e < 1 — ctie, 
for sufficiently small e with ct\ = y/ ei/(2E ) and a 2 = 2E /e 2 . □ 

5. Numerical Results 

In this section, we will present the results of numerical simulations that illustrate and 
clarify our analysis of the energy and also the structure of the minimizers for the regularized 
Cross-Newell energy jF e within the class of functions given by (2). 

For our numerical simulations, we restrict ourself to the finite domain, lZ e = {(x, y) \ < 
x < I = 7r/e,0 < y < L}, where L 3> 1 is a length scale much larger than the typical 
wavelength of the pattern. The boundary conditions in (2) which are appropriate for the 
semi-infinite strip S € are modified for the finite domain as follows - 

9(x, 0) = for < x < al; 

y (x, 0) = for al < x < I; 

9(x, y) — ex is periodic in x with period I for each y e [0, L]; 

(47) 9(x, L) = \ex + y/l - e 2 L + 5 

It is rather straightforward to show that there exist 9 e e H 2 (1Z e ) satisfying (47) for an 
a e G [0, 1) and 5 e G R minimizing the functional 

jF e [0;a,5] = J J {[A9] 2 + (l-\V9\ 2 ) 2 }dxdy. 
11 22 



The existence of a minimizer is immediate from the following lemma: 



Lemma 5.1. Let < a < 1, and 5 G K be given, pj G L 2 (1Z e ) is a sequence of functions 
that converges weakly to zero. iJp er denotes the completion of periodic (in x) functions on 
1Z e with respect to the H 2 norm. If 9j G H 2 er (1Z e ) is a sequence satisfying (in the sense of 
trace) 

A9j = Pj 
6 j (x,0) = forO<x <al; 
d y 8j(x, 0) = for al < x < I; 
(48) 9(x,L)=0 

it follows that, up to extraction of a subsequence and relabelling, we have V9j — > in 

L 4 {n e ,R 2 ). 

Proof. Elliptic regularity along with the given boundary conditions implies that the sequence 
9j is bounded in H 2 (1Z e ). The compactness of the embedding H 2 (1Z e ) <^-> W 1,4 (1Z e ) [7] proves 
the lemma. □ 

If 9j is an infimizing sequence for T e \9]a, 5} subject to the boundary conditions in (47), 
then let 9j = 9j — (p, where ip is a smooth function on lZ e satisfying the boundary conditions 
in (47). It then follows from the form of jF e and the fact that 9j is infimizing that A9j is a 
bounded sequence in L 2 , and so converges weakly to a limit p*. Applying the compactness 
result of the preceding lemma with reference to the sequence pj = A9j — p*, we obtain the 
existence of a minimizier for the functional ^[9; a, 5] for a fixed a and 5. 

Note that, for a given a it is easy to construct smooth transformations ip t '■ 7t e — y 7Z e such 
that ipQ is the identity, if 9 satisfies the boundary conditions in (47), then 9 o ip t satisfies 
the same boundary conditions with the fraction of the boundary with a Dirichlet boundary 
condition equaling a(l + t). Further, the energy F e [9 o tp t - a(l + t), 5] is a smooth function 
of t for sufficiently small t. A standard argument now implies that, for a given S the map 

a i— > inf JF e [#; a, 5} 

is continuous for a G (0, 1). A similar argument shows that the map is also continuous at 
a = 0. In Lemma 4.1 we showed that liminf a ^! T^\9\ a, 5] = oo. Combining these results, 
we see that 

miF e [9;a,5] = inf 

a,8 a 

We now consider variations 9 — > 9 t = 9 + tx{y/L), where x is a smooth function vanishing 
identically on [0,1/3] and equal to 1 on [2/3,1]. The functions 9 t satisfy the boundary 
conditions in (47), except the asymptotic phase shift is given by 5 + t. A similar argument 
as above shows that the map 

5 ^ miF £ [9;a,5] 

is continuous and it is easy to see that T^\9\ a, 5] — > oo as 5 — > ±oo. In particular, this 
proves the existence of an optimal 5, and combining with the results from above, we see that 
the minimizer 9 e , a e , S e , can be obtained by successive minimization in each of the factors. 
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inf F £ [9;a,5} . 



This suggests the following discretization for the functional T e , which should converge as 
the grid spacings 77, ^" — ^ 0. We define a grid by Xi = irj,i = 0,1,2, ... ,m — 1, yj = = 
0,1,2, ... ,n, where i] = l/m,( = L/n. The discretization of the test function 6{x, y) is 

We define the difference operator 5^ by 

&) M = ± w- g * 

with similar definitions for 5^. In terms of the discretization, the boundary conditions are 

iiO = for < % < k; 
3y9i,o = for k < i < m; 

(49) nhj =6 0tj + 7t j = 0,1,2,..., n 

(50) 9i n = — + y/l - e 2 L + 5 i = 0, 1, 2, . . . m - 1 

m 

and the Energy functional is discretized as 

"(M + + (W?, + (5-9)1 



m— 1 n 



i=0 j=0 



- 1 



Computing this functional requires assigning values for 6ij with i = —l,j = —1 and j = n+1. 
The values for i — — 1 are obtained form the shift-periodicity of 9 by 9-ij = 9 m -\j — ir. 
The values for T j —n + 1 are assigned using the Dirichlet boundary condition 9 iyn+ i = 
^ + y/l — e 2 (L + + 5. This functional is minimized using MATLAB's conjugate-gradient 
minimization. 

Fig. 3 shows the results from minimizing the RCN energy over the pattern 9 and also 
the phase shift 5, for different values of e, and for a range of values of a k/m. The 
results do indeed suggest that the (partial) minimization with respect to the pattern and the 
asymptotic phase yields a functional that depends continuously on a. Further, this functional 
has first-order (discontinuous) phase transition at a bifurcation value e*, below which the 
global minimizer has a 7^ 0. 

Fig. 4 shows the energy of the minimizer (minimizing over the pattern, asymptotic phase 
and the parameter a) as a function of e. Note that the minimum energy is a non-differentiable 
function of e, as one would expect for a first-order phase transition. 

Figure 5 show the numerically obtained minimizing patterns at various values of e. Note 
that, for sufficiently large e, the minimizers are the knee-solutions (3) with a = 0, whereas for 
sufficiently small e, the minimizers have convex-concave disclination pairs, and have a ^ 0. 
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Figure 4. Global minimum of the regularized Cross- Newell energy. 
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